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^ ■ ABSTRACT 

^ ' Aims. We determine the effective temperature, surface gravity and projected rotational velocity of 

nine T dwarfs from the comparison of high-resolution near-infrared spectra and synthetic models, 
and estimate the mass and age of the objects from state-of-the-art models. 
Methods. We use the AMES-COND cloudless solar metallicity models provided by the 
. PHOENIX code to match the spectra of nine T-type field dwarfs observed with the near-infrared 

high-resolution spectrograph NIRSPEC using ten echelle orders to cover part of the J band from 
1.165 to 1.323 yum with a resolving power R -20,000. The projected rotational velocity, effec- 
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^\ ^ tive temperature and surface gravity of the objects are determined based on the minimum root 

mean square of the differences between the modelled and observed relative fluxes. Estimates of 
the mass and age of the objects are obtained from effective temperature- surface gravity diagrams, 
' where our results are compared with existing solar metallicity models. 

■ Results. The modelled spectra reproduce quite well the observed features for most of the T 

dwarfs, with effective temperatures in the range of 922-1009 K, and surface gravities between 
lO"*"^ and 10^" cms^-. Our results support the assumption of a dust free atmosphere for T dwarfs 
later than T5, where dust grains form and then gravitationally sediment into the low atmosphere. 
The modelled spectra do not accurately mimic some individual very strong lines like the K i dou- 
blet at 1.2436 and 1.2525 ^m. Our modelled spectra does not match well the observed spectra 
of the two T dwarfs with earlier spectral types, namely SDSSp J125453.90-012247.4 (T2) and 
2MASS J05591914-1404488 (T4.5), which is likely due to the presence of condensate clouds 
that are not incorporated in the models used here. By comparing our results and their uncertain- 
ties to evolutionary models, we estimate masses in the interval s;5-75 Mj for T dwarfs later than 
T5, which are in good agreement with those found in the literature. We found apparent young 
ages that are typically between 0.1 and a few Gyr for the same T dwarfs, which is consistent with 
recent kinematical studies. 



Key words, atomic lines - molecular lines - low-mass objects - ultracool dwarfs spectra 



2 Please give a shorter version with: \authorruiming and/or \titilerunning prior to \maketitle 

1. Introduction 

The first near-infrared spectrum of a T dwarf (Gl 229B) showed a predominance of methane (CH4) 
bands that made it look similar to that of Jupiter (Oppenheimer et al. 1995) despite of a difference 
in temperature of about 800 K. Atmospheric models were soon developed to fit the spectrum of 
Gl 229B to estimate its surface gravity, effective temperature, age and mass (Allard et al. 1996; 
Marley et al. 1996; Tsuji et al. 1996). These models indicated that the atmosphere of this brown 
dwarf was free from dust grains, which were needed to explain the properties of the warmer L 
dwarfs (Allard et al. 1997). The dust grains in Gl 229B could be condensed in clouds dominated 
by organic compounds (Fegley et al. 1996; Griffith et al. 1998). It has also been proposed that a 
warm dust layer could be present deep inside the photosphere (Tsuji et al. 1999) and that the alkali 
resonance lines have very strong pressure-broadened red wings that provide a significant source of 
opacity at near-infrared wavelengths (Burrows et al. 2000). 

The discovery of free-floating objects similar to Gl 229B (Burgasser et al. 1999; Cuby et al. 
1999; Strauss et al. 1999) indicated that these objects are numerous in the solar vicinity. A uni- 
fied near-infrared classification scheme for T dwarfs has been defined by Burgasser et al. (2006). 
Theoretical models provide a good correspondence to the broad-band colors and low-resolution 
spectra of T dwarfs (Burrows et al. 2006). Currently, the coolest T dwarfs known have estimated 
effective temperatures between 600 K and 700 K (Warren et al. 2007; Delorme et al. 2008). 

Most modelling efforts on T dwarfs have concentrated on fitting broad-band colors and low- 
resolution spectra, and deriving atmospheric parameters by comparing the data with synthetic spec- 
tra (Gushing et al. 2008, Leggett et al. 2007, Saumon et al. 2007, Saumon et al. 2006, Burgasser 
et al. 2006, Tsuji et al. 2005, Tsuji et al. 2004, Burgasser et al. 2004). Martm & Zapatero Osorio 
(2003) estimated the surface gravity and effective temperature of one T dwarf from a mid-resolution 
near-infrared spectrum. High-resolution (R~20,000) observations in the /-band obtained with 
NIRSPEC on the Keck II telescope are now available (Zapatero Osorio et al. 2006, Mc Lean et 
al. 2007). 

Modelling the atmospheres and spectra of T dwarfs presents a set of challenges. In the AMES- 
COND models used in this work (Allard et al. 2001), the profile of individual strong IR lines are 
approximated by Voigt profiles with estimated damping constants as no better data are currently 
available for these lines. The line databases for molecules such as FeH, CaH, CH4 are not complete 
and in many cases lack accuracy. However, in general the molecular line data are quite good as data 
sources such as Hi l RAN are of high quahty. The lower the effective temperatures are, the better the 
completeness of the molecular fine data as the (comparatively) less accurate higher energy levels 
are less populated and thus have smaller effects at lower temperatures. 

In the AMES-COND models a dust-free configuration is used, where the dust particles form 
but do not contribute to the opacities (Allard et al. 2001). This assumption should be reasonable 
for cool T dwarfs, but will become progressively worse for effective temperatures above 5»1200 K. 
Overall, the systematic errors are still considerably larger than, e.g., for solar type stars, but the 
situation can only improve if more comparison between models and data are made to identify the 
areas where improvements are most urgently needed. 
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In this paper, we present the first comparison with AMES-COND models of high-resolution 
near-infrared spectra of a sample of nine T dwarfs observed by Zapatero Osorio (2006) that are 
used to derive their physical parameters. A summary of the observations used here is given in 
Section |2] Section [3] is devoted to the analysis of the data. The results are presented in Section |4l 
Section|5]contains the discussion of our results, and finally. Section |6]contains our conclusions. 

2. Observations, models and data preparation 

2.1. Observations 

High-resolution near-infrared spectra in the 7-band of nine T dwarfs were obtained using the 
Keck II telescope and the NIRSPEC spectrograph. The instrumental setup was chosen to provide a 
wavelength coverage from 1 . 148 up to 1 .346 fim spht into ten different orders, a nominal dispersion 
ranging from 0.164 (blue wavelengths) to 0.191 Apix"' (red wavelengths), and a final resolution 
element of 0.55-0.70 A at 1.2485 fim (roughly the central wavelength of the spectra), correspond- 
ing to a resolving power R - \1 , 800 - 22, 700. A detailed description of the observations and data 
reduction is provided in Zapatero Osorio et al. (2006) (hereinafter ZO06), who presented the T 
dwarfs that we used in this study. 

In our spectroscopic analysis we have used as many echelle orders as possible: from order 
57 up to order 66; this depends on the quality (in terms of signal-to-noise ratio) of the observed 
spectra. Table 2 in McLean et al. (2007) summarizes the wavelength and dispersion properties 
of each of the /-band NIRSPEC echelle orders 58 to 65 for the same instrumental configuration 
that we have used here. Order 57 covers the wavelength interval 1.327-1. 347 jum with a nominal 
dispersion of 0.191 Apix Order 66 covers the range 1.147-1. 164 yum with a nominal dispersion 
ofO.167 Apix-i. 

Fig. [U shows, as an example, the order 61 with the presence of the Ki doublet at 1.2436 and 
1.2525 //m, for all the T dwarfs. 

2.2. Description of atmospliere modeis 

The PHOENIX code (Hauschildt & Baron 1999) is a general purpose stellar atmosphere modelling 
package that includes very complex atomic models and line blanketing by hundreds of millions 
of atomic and molecular hnes. The code can be used to compute model atmospheres and synthetic 
spectra of cool objects. The radiative transfer in PHOENIX is solved in spherical geometry using an 
operator splitting method (Hauschildt 1992, 1993). At the low temperatures of the brown dwarfs, a 
rich chemistry in the atmospheres of these objects is present, with hundreds of gas-phase species, 
liquids and crystals, and the formation of tens of different types of dust grains (e.g., silicates, 
amorphous carbon, iron). There are four different scenarios for the dust formation considered in 
PHOENIX, among which is the AMES-COND cloudless models (AUard et al. 2001). 

Here we used a grid of the AMES-COND v2.2 models and synthetic spectra of solar metallicity, 
with effective temperatures Jeff ranging from 700 to 3000 K (steps of lOOK) and surface gravities 
log g ranging from 3.0 to 5.5 (steps of 0.5) with g in cm s^^, to mimic the observed spectra of T 
dwarfs. 
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SDSSp J125453.90-012247.4 T2 



1.245 1.250 1.255 

wavelength (jum) 

Fig. 1. Variation of the strong K i doublet at 1 .2436 and 1 .2525 fim along the spectral type sequence 
of T dwarfs. All spectra are normalized to unity over the entire wavelength range and offset by a 
constant. Also, each spectrum is shifted in velocity to vacuum wavelengths. 



The AMES-COND v2.2 models use the H2O line list of Barber et al. (2006), the FeH list of 
Dulick et al. (2003), and the damping constants used for the K i lines are published in Allard et al. 
(2003). 

Fig.|2]illustrates the identification of some near-infrared features over a synthetic spectrum with 
Vrot sin/ - 0, Tetf = 1000 K and log g = 4.5. Most features are due to water vapor, and the only visible 
atomic lines are due to the K i doublet at 1 .25 ;um. 

Fig. [3] shows the variations of the synthetic spectra for a surface gravity of log g - 4.5 and 
dififerent effective temperatures from 700 to HOOK. We note the strong change of the Ki doublet 
at 1.2436 and 1.2525 fim with temperature. Fig.|4]shows the variations of the synthetic spectra for 
various log g (from 3.5 to 5.5) and the same Teif^lOOOK. 

2.3. Data preparation 

All of our observed spectra were moved to vacuum wavelengths (i.e., laboratory frame of refer- 
ence) for a proper comparison with the theoretical models. This was done from a cross-correlation 
analysis. In a first iteration, we used the radial velocities published in Zapatero Osorio et al. (2007). 
Then, we determined the atmospheric parameters for each T dwarf as explained later and used the 
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Fig. 2. Synthetic flux versus wavelength. Identifications of some features of a synthetic spectrum 
with Teff =1000 K and log g = 4.5. Most of the Unes are due to water vapor. 



corresponding synthetic spectra as templates to recompute the radial velocities. New atmospheric 
parameters were thus obtained consistently with the recalculated radial velocities. We note that our 
radial velocities are in agreement, within the errors, with those of Zapatero Osorio et al. (2007). 

The grid of PHOENIX synthetic spectra were modified in order to be compared with the 
NIRSPEC observations. First, the synthetic spectra were transformed to take into account the pro- 
jected rotational velocity (Vrot sinj) of the objects using the formaUsm of Gray (1992), with a limb 
darkening parameter e=0.6. A grid of models with projected rotational velocities between 10 and 
50 Km s"\ with steps of 1 Km s~\ were generated. These spectra were also convolved with a 
Gaussian that mimics the instrumental profile along the dispersion axis for each order. The result- 
ing modelled spectra were rebinned to the same resolution of the observations. Modelled spectra are 
normalized over the wavelength range corresponding to order 59. All these steps were performed 
using our own programmes. 
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Fig. 3. Synthetic flux versus wavelength. Models with logg = 4.5 (CGS units) and T^b^IOO, 800, 
900, 1000 and 1 100 K. The synthetic spectra were computed with 0.05 Apix~\ but are shown after 
being smoothed with a boxcar average of 10 pixel width. 
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Fig. 4, Synthetic flux versus wavelengfli. Models with reif=1000K and logg = 3.5, 4.0, 4.5, 5.0 
and 5.5 (CGS units). The synthetic spectra were computed with 0.05 Apix~\ but are shown after 
being smoothed with a boxcar average of 10 pixel width. 

3. Analysis 

In order to constrain the number of possible solutions provided by our large set of models, the 
root-mean- square RMS (Vrotsim, T^g, log g) is obtained for each model: 



RMS (Vrotsini, T^eff , log g) 



i 



(1) 
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Table 1. SDSSp J125453.90-012247.4. Synthetic and observed spectra comparison for all orders: 
RMS, N, Teff, log g and Vi-ot sin/. 



order N RMS T^g logg sini 



00 


yyo 


u.z/y / 


loUU 


c c 
J.J 


1 O 
JO 


OJ 


yyj 


U. i /4U 


1 ^nn 


c c 

J.J 


40 


O'^ 


1 nni 


U. 1 jJiO 


9 Ann 


J.J 


jO 


63 


1002 


0.0931 


2500 


3.5 


26 


62 


1002 


0.0839 


2700 


3.5 


27 


61 


1002 


0.0960 


2700 


5.5 


46 


60 


890 


0.0914 


700 


5.5 


28 


59 


1002 


0.0811 


2000 


3.0 


36 


58 


1002 


0.0932 


1500 


5.5 


36 


57 


994 


0.0862 


2000 


4.0 


31 



Table 2. 2MASS J05591914-1404488. Synthetic and observed spectra comparison for all orders: 
RMS, N, Teff, logg and Vrot sin/. 



order 


N 


RMS 




logg 


v,ot sini 


66 


987 


0.1079 


1500 


5.0 


21 


65 


989 


0.0796 


1200 


5.5 


25 


64 


991 


0.0575 


1000 


3.5 


33 


63 


991 


0.0512 


800 


5.5 


24 


62 


991 


0.0416 


800 


5.5 


27 


61 


991 


0.0612 


1000 


5.0 


33 


60 


882 


0.0434 


800 


5.5 


25 


59 


991 


0.0724 


900 


3.5 


33 


58 


991 


0.0682 


1300 


5.5 


33 


57 


957 


0.0491 


1500 


4.0 


33 



where / stands for i-pixel in the spectral axis, Fy(i) and Gv(i) correspond to the observed and 
modelled fluxes, respectively. We compute values of RMS for the whole set of models and for all 
echelle orders independently, since different orders show diff'erent signal-to-noise ratios. 

Tables [T]|9] show the values of T^ff, log g, and v'rot corresponding to the model that minimizes 
the RMS for each echelle order. The total number of pixels used in the comparison is also listed. 
For order 60, the wavelength range between 1.268 and 1.270 ;um is excluded from the comparison 
because of the presence of strong telluric lines that practically blocked any signal from the targets. 
To avoid border effects, a few pixels at both sides of each echelle order have been rejected. 

3. 1 . Average physical parameters 

We compute the average values of Teff, log^, and Vrot for each object in our sample by combining 
the various values obtained for the different echelle orders according to the following equations: 
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Table 3. 2MASS J1503 1961 +2525 196. Synthetic and observed spectra comparison for all orders: 
RMS, N, Teff, log g and v,-ot sin/. 



order 


N 


RMS 


Tff 


los 3 
^'^o o 


Vivit sin/ 


66 


969 


0.1112 


1300 


5.0 


34 


65 


978 


0.0627 


1100 


4.0 


42 


64 


987 


0.0542 


900 


4.0 


30 


63 


987 


0.0445 


1000 


3.5 


42 


62 


987 


0.0413 


800 


5.5 


37 


61 


987 


0.0576 


1000 


5.5 


42 


60 


873 


0.0461 


800 


5.5 


31 


59 


987 


0.0619 


900 


3.5 


37 


58 


987 


0.0633 


1300 


5.0 


33 


57 


911 


0.0448 


1400 


4.5 


33 



Table 4. SDSS J162414.37+002915.6. Synthetic and observed spectra comparison for all orders: 
RMS, N, Teff, logg and Vrot sin/. 



order 


N 


RMS 




logs 


Vi-o, sin; 


64 


977 


0.0396 


1000 


4.0 


43 


63 


977 


0.0437 


1000 


4.5 


47 


62 


977 


0.0386 


800 


5.0 


41 


61 


976 


0.0550 


900 


5.5 


47 


60 


858 


0.0533 


900 


3.5 


37 


59 


977 


0.0582 


1100 


5.0 


40 


58 


977 


0.0671 


1300 


5.0 


35 



Table 5. SDSS J134646.45-003150.4. Synthetic and observed spectra comparison for all orders: 
RMS, N, Teff, logg and Vrot sin/. 



order 


N 


RMS 


Teff 


logg 


Vrot sini 


64 


972 


0.0544 


1000 


4.0 


19 


63 


974 


0.0554 


1000 


4.0 


20 


62 


974 


0.0500 


900 


4.5 


17 


61 


974 


0.0644 


1000 


5.0 


25 


60 


862 


0.0801 


900 


4.0 


12 


59 


974 


0.1185 


900 


4.0 


13 


58 


973 


0.1139 


1200 


4.5 


13 



<logg>^ ly. ■ (3) 

a-j 

2y(v,-otsin/;^) 

< Vi-otsin/ >= , (4) 

V — i 

a-j 

where crj refers to the minimum RMS for the echelle order j (as hsted in Tables [T]|9]l, and Wj 
corresponds to the weight that we have assigned to the echelle order j based on the signal-to-noise 
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Table 6. 2MASS J15530228+1532369. Synthetic and observed spectra comparison for all orders: 
RMS, N, Teff, logg and Vi-ot sin/. 



order 


N 


RMS 




log g 


Viot sini 


63 


976 


0.0507 


1000 


4.5 


37 


62 


976 


0.0582 


800 


5.0 


21 


61 


976 


0.0659 


900 


5.5 


39 


60 


857 


0.0563 


900 


4.0 


27 


59 


976 


0.0733 


900 


4.0 


30 


58 


979 


0.0863 


1100 


5.0 


23 



Table 7. 2MASS J12171110-0311131. Synthetic and observed spectra comparison for all orders: 
RMS, N, Teff, log g and v,-ot sin/ 



order 


N 


RMS 




logs 


v,„t sini 


63 


843 


0.0467 


800 


5.5 


21 


62 


972 


0.0581 


900 


4.0 


32 


61 


972 


0.0947 


900 


5.5 


41 


60 


882 


0.1777 


900 


4.0 


21 


59 


972 


0.0894 


1000 


4.5 


27 


58 


972 


0.0918 


1100 


4.5 


29 



Table 8. GL 570D. Synthetic and observed spectra comparison for all orders: RMS, N, Tes, log g 
and v'rot sin/. 



order 


N 


RMS 


Tts 


logg 


Vi-ot sini 


64 


974 


0.0327 


900 


4.5 


27 


63 


974 


0.0411 


1000 


4.5 


34 


62 


974 


0.0428 


1000 


4.0 


41 


61 


974 


0.0649 


900 


5.5 


41 


60 


861 


0.0675 


900 


4.0 


26 


59 


974 


0.0919 


1000 


4.5 


26 


58 


974 


0.0986 


1000 


5.0 


21 



ratio of the observations. TablefTOlshovys the various Wj, which are relative to order 61. Regarding 
the echelle order 60, we have given lower weights to J 12 17-03 11 and J0415-0935 because the 
correction for Earth telluric lines were not optimal. Finally, the average values of T^ff, log g, and 
Vrot are provided in Table (TT] where we also include the spectral types of the targets following the 
unified scheme of Burgasser et al. (2006). Errors in < T^ff >, < log g > and < Vjotsin/ > were 
obtained from the standard deviation of the corresponding values derived for the various echelle 
orders. We note that if no weights are taken into account, the solutions for < T^ff >, < log g > and 
< v'rotsin/ > do not change significantly, with diff'erences with respect of those obtained using the 
weights Wj of +40 K (for T dwarfs later than T5 the difference is only ±7 K), +0.1 dex, and +1 
Km s ^ which are much smaller than the quoted uncertainties. 

Fig. |5] shows, as an example, some contours of equal RMS (Root mean square) around the 
solution (< Teff >, log g) for GL570D, with v,ot=32 Km s"'. The isocontour of the RMS of 1.25 
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Table 9. 2MASS J0415 1954-0935066. Synthetic and observed spectra comparison for all orders: 
RMS, N, Teff, log g and Vi-ot sin/. 



order 


N 


RMS 


Tff 


■^^6 o 


Vpot sin/ 


64 


926 


0.2029 


800 


3.5 


43 


63 


981 


0.1791 


1000 


3.5 


33 


62 


982 


0.1119 


1000 


4.0 


44 


61 


987 


0.1206 


900 


5.0 


43 


60 


899 


0.1539 


1000 


3.5 


38 


59 


987 


0.1081 


900 


4.5 


29 


58 


984 


0.1330 


1000 


5.0 


25 



Table 10. Weights of echelle orders. 



Order 


Weight 


Remark 


66 


0.3 


for all objects 


65 


04 


for all objects 


64 


0.6 


for all objects 


63 


1.0 


for all objects 


62 


1.0 


for all objects 


61 


1.0 


for all objects 


60 


0.2 


for2MASS J12171 1 10-031 1 131 




0.6 


for 2MASS J0415 1954-0935066 




1.0 


for the rest 


59 


1.3 


for all objects 


58 


1.1 


for all objects 


57 


0.8 


for all objects 



times the minimum RMS enclose values of logg of ^ +0.7 dex and Teff of ~ +150K around the 
solution. In general, the isocontours of the RMS enclose a unique solution with a broad range of 
values of gravity and temperature around the solution. 

3.2. Best modelled spectra 

We have computed the modelled spectra determined from the average effective physical parameters 
determined in Sect. 13 . II (hereinafter, best modelled spectra). We have used a bilinear interpolation 
to compute the synthetic models coiTesponding to < Teff > and < log g > from om AMES-COND 
model grid. This bilinear interpolation has been done using the logarithm of the synthetic fluxes. 
We first proved that a linear interpolation in both axes was appropriate. We fixed T^ff at lOOOK and 
used the corresponding synthetic models at log ^=4.0 and log g=5 .0 to compute interpolated spectra 
at logg=4.5. We also fixed logg=4.5 and used the coiTesponding synthetic models at Teff =900K 
and HOOK to compute the interpolated spectra at T^ff =1000K. We then applied the same steps 
described in Sect 2.3 to compare the synthetic spectra to our observations. The diff'erences between 
the interpolated and synthetic models are typically of a few per cent. The largest differences of 
10% observed between the spectra were obtained for the echelle order 57 with a null rotational 
broadening. For the rest of the echelle orders the differences were generally below 1%. 



12Please give a shorter version with: \authorrunning and/or \titilerurming prior to \maketitle 




3.0 1 . I ... I ... I ... I ... I ... I .. . 

800 1000 1200 1400 1600 1800 2000 
T [K] 

Fig. 5. Contours with RMS=1.25 (thick Hnes) and 1.5 (thin lines) times the value of the minimum 
RMS for orders 59 (solid lines) and 60 (dashed lines). Location of (< Teff >, log g) for the minimum 
RMS is also indicated (dot). 



Table 11. Projected rotational velocity, temperature, surface gravity, mass and age of T dwarfs. 



Object 


Spectral Type 


Vrotsini 




PREV. 


REF. 


< logg > 


PREV. 


REE 


M 


Age 








K 


K 




(cm s^^) 


(cm s^^) 




Mj 


Gyr 


SDSSp J125453.90-012247.4 


T2 


34±7 


2007+662 


1200 


1 


4.3+1.0 


5.0 


1 






2MASS J05591914-1404488 


T4.5 


28±5 


1002+278 


1200 


1 


4.9+0.8 


5.5 


1 


30:« 


1 0+' 


2MASS J 1503 196 1+2525 196 


T5 


36±5 


1009+217 


1178 


2 


4.6+0.8 


«5.0 


5 






SDSS J162414.37+002915.6 


T6 


42±5 
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Figs.l6]fl4lshow the observations and modelled spectra for our list of T dwarfs. 
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4. Results 

4.1. Near-infrared absorption features of T dwarfs 

As observed by McLean et al. (2007), the 7-band spectral morphology at R ~ 20, 000 of T dwarfs 
shows a dense population of weak absorption features and a few relatively strong lines. The fine- 
scale spectral structure is mainly associated with H2O, and the strong atomic features are doublets 
of Ki that appear in orders 61 and 65. The strong Ki doublet in order 61, at wavelengths 1.2436 
and 1.2525 //m, deepen towards earlier Ts (see Fig. [T]i. We note that the lines of this doublet in 
SDSS J 134646.45-003 150.4 are stronger than what would be expected according to its spectral 
type (T6.5). 

FeH hnes in orders 62 and 63 are weaker towards later Ts. Some absorption around 1.222 fim 
may be present even at T4.5. Those at 1.24637 and 1.24825 i-im FeH in order 61 are not visible in 
T dwarfs. Order 64 is dominated by sharp and deep H2O absorption features. This is also the case 
of order 66, and more dramatically for order 57, with the strongest H2O absorption in the observed 
range. Order 65 shows the other Ki doublet and changes with spectral type. It also contains many 
strong intrinsic transitions of hot H2O, like the feature at 1.175 fim. 

In general, our observed spectra are well reproduced by the best modelled spectra. In a few 
cases we have found some small differences. This is the case of 2MASS J15031961H-2525196 (T5) 
(see Fig. |8j, where the best model is found for log g-4-.6 and T^ff- 1 009 K. The observed K i doublet 
in order 61 is not perfectly reproduced by the model. We also note some differences in orders 59 
and 58. Small differences in these orders are also present in SDSS J162414.37H-002915.6 (T6) (see 
Fig.©. These are much smaller in later spectral types (e.g., SDSS J134646.45-003150.4 (T6.5), 
see Fig. [Toll. For late T dwarfs, with faint K i doublets, the models provide an excellent match to the 
observations for all the orders, with a few exceptions such as the order 61 in 2MASS J12171110- 
0311131 (see Fig.[T2li. where the model is flat towards the blue region. We note that the observed 
spectrum may be affected by errors in the flat-fielding. Note also that there is a faint line of Ki at 
1.2525 yum in the best modelled spectra for order 61 of GL 570 D that disappears when considering 
an effective temperature 48 K lower (see Fig.flJli. 

The two earlier T dwarfs of our sample, SDSSp J 12545 3. 90-0 12247 .4 and 2MASS J05591914- 
1404488, present significant differences between the best modelled and observed spectra. The 
modelled spectrum corresponding to the average values (reff=1002K and log g -4.9) in 2MASS 
J0559 19 14- 1404488 fails to mimic order 57 with many water vapor lines. The Ki doublet at 1.2436 
and 1 .2525 //m, very sensitive to temperature, is also difficult to reproduce. Order 57 is clearly bet- 
ter modelled with a higher temperature Teff-llOOK. 

4.2. Pliysical parameters of T dwarfs 

4.2.1. Rotation, effective temperature and surface gravity 

For our list of nine T dwarfs, we have determined their projected rotational velocity (vrotsini), ef- 
fective temperature and surface gravity from the comparison of the high resolution spectra with 
the AMES-COND models (see Table fTTT i. The rotational velocities determined here from mod- 
elled spectra are in good agreement (i.e., within l-cr uncertainties) with those obtained by ZO06, 
who used J134646. 45-003150.4 as a non-rotator reference to measure the Vrotsin/ parameter of the 
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sample (i.e., they assumed Vi-otsin/=0 for this object). With the exception of 2MASS 12171110- 
031 1 131, our values are, indeed, systematically higher by a few Km s which is expected since 
the template used by ZO06 has a small rotation of 16 Km s ' according to our analysis. 

The effective temperatures of the T dwarfs with spectral types later than T5 are between 922 
and 1009 K, with errors between ^50 and 200 K (see Table fTTTi. The surface gravities of those 
objects lie between 10"* ' and 10'^-^ cms^^, and errors are of aiO.7 dex. In Sections |5 .21 and |53] we 
discuss our results on Teff and log g. 

4.2.2. Mass and age 

We have estimated the mass and age of our T dwarfs using state-of-the-art evolutionary models. Our 
spectroscopically derived values < Teff > and < logg > are plotted in the two panels of Fig. [15] The 
T2V dwarf J 1254-0 122 is discarded from the Figure since the modelled spectra found here do not 
reproduce the observations satisfactorily. Therefore, only mid- and late-T dwarfs are shown in the 
Jeff versus logg diagram, which we remark is independent of the distance to the sources. Fig. [15] 
{left) shows, overplotted on the data, the solar metallicity models by the Lyon group (Baraffe et 
al. 2003), with isochrones (from lOMyr to lOGyr) and curves of constant mass (5, 10, 20, 30, 
and 50 Mj). Fig.[T5](Wg/if) shows our data with the solar metallicity models by the Arizona group 
(Burrows et al. 1997): isochrones (10 Myr, 0.1, 1 and lOGyr) and curves of constant mass (5, 10, 
20, 30, 50, 70 and 75 Mj). 

Table [TT] shows the values of mass and age for our sample of T dwarfs that are derived from the 
Teff versus log g diagram using the models of the Lyon group. The mass and age uncertainty inter- 
vals are inferred from the error bars in the effective temperature and surface gravity. Summarizing, 
all these objects have a mass in the interval 5-75 Mj, thus confirming their very likely substellarity, 
and likely ages that interestingly seem to be younger than the solar system. We note, however, that 
considering the large error bars, most of the T dwarfs in our sample may have an age consistent 
with a few to several Gyr. For GL570D and 2MASS J0415 1954-0935066 we have derived an age 
upper limit of 1-2 Gyr, and only for SDSS J 134646.45-003 150.4 the upper Umit is below 1 Gyr. 

5. Discussion 

5. 1 . Modelled and observed spectra 

Most of the lines in a T dwarf are due to water vapor. A few atomic lines can be identified, but 
some of them, e.g. the Rb I line, are very weak compared to the haze of molecular lines. 

In general, the modelled spectra match remarkably well the high-resolution observed spectra. 
We note only that the profiles of the strong K i lines are not so well matched by the models. 
This is likely due to the approximations used for the intrinsic line profiles in the models, most 
importantly estimates of the damping constants. For example, Allard et al. (2003) and Johnas et 
al. (2007) have demonstrated the importance of detailed line profiles (not just damping constants) 
for optical alkali resonance lines in cool sub-stellar objects, however, such profiles are not yet 
available for the IR alkali lines. Therefore, the profiles of these atomic lines cannot be modelled 
accurately, and so cannot be used to obtain reliable parameter estimates. From qualitative con- 
siderations, however, the models predict that for surface temperatures of 900-lOOOK, the atomic 
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Fig, 6. SDSSp J125453.90-012247.4. Black and red lines correspond to the observed and modelled 
spectra, respectively. The model is for an effective temperature of 2000 K and a surface gravity of 
lO"*-^ cm s-2. 



lines of Ki strengthen with decreasing gravity, as was pointed out by Knapp et al. (2004). SDSS 
J134646.45-003150.4shows the Ki doublet at 1.25 yum, in contrast to SDSS J162414.37+002915.6 
and 2MASS J15530228+1532369, which have similar spectral types. From our analysis, we have 
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Fig. 7. 2MASS J05591914-1404488. Black line corresponds to the observed spectrum, and red and 
blue lines correspond to modelled spectra with the same g=lO'^-^ cm s"^ and reff=1002, and 1700 
K. 



found the smallest log g for SDSS Jl 34646.45-003 150.4 in our sample. This is also the object with 
the lowest estimated age. 

The FeH lines are difficult to model, in part due to their relatively poorly known oscillator 
strength and possible errors in the chemistry. The fact that they become weaker with lower ef- 
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Fig, 8. 2MASS J1503 1961 +2525 196. Black and red lines correspond to the observed and modelled 
spectra, respectively. The model is for an effective temperature reif=1009 K and g-lQ'^-^ cm s"^. 



fective temperatures (like those of our sample) is due to the condensation of iron, which reduces 
the number density of FeH, resulting in an important test of the treatment of condensation in the 
equation of state. 
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Fig. 9. SDSS J162414.37+002915.6. Black and red lines correspond to the observed and modelled 
spectra, respectively. The model is for an effective temperature of 980 K and a surface gravity of 
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CH4 bands are expected and have been observed in T dwarfs (Burgasser et al. 2002) but, in 
the spectral range considered here, are much weaker than the water vapor lines and so cannot be 
identified with confidence. 

At higher effective temperatures, as for SDSSp J125453.90-012247.4 (T2) and 2MASS 
J05591914-1404488 (T4.5), the models appear to be generally less consistent with the observa- 
tions. This is likely due to the presence of remnants of dust clouds, floating around in upper layers, 
that are not included in the modelled spectra (see Ruiz et al. 1997, Ackerman & Marley 2001, 
Burgasser et al. 2002, Burrows et al. 2006, Cooper et al. 2003, HeUing et al. 2008). Gushing et al. 
(2008) have computed the properties of SDSSp J125453.90-012247.4 and 2MASS J05591914- 
1404488 from the comparison of low and intermediate resolution spectra in the 0.95-14.5 fim 
wavelength range and synthetic spectra. They found that SDSSp J125453. 90-012247. 4 has con- 
densate clouds that are thicker than those in 2MASS J05591914-1404488, which may explain the 
difference in the spectra of these overluminous T dwarfs. Gushing et al. also discussed the possible 
unresolved binarity in both objects. They determined Teff ~1200K (<1150K), log g ~4.8 (<5.38) 
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Fig. 10. SDSS J 134646.45-003 150.4. Black and red lines correspond to the observed and modelled 
spectra, respectively. The model is for an effective temperature of 990 K and a surface gravity of 
cm s" 



and an age of -0.4 (<10) Gyr for 2MASS J05591914-1404488 assuming it is a single object (equal 
mass binary). The same temperature of ~1200K is found for SDSSp J125453.90-012247.4. Our 
values of Teff, log g and age of 2MASS J05591914-1404488 (see Table[Il]i are in agreement with 
those of Gushing et al. 

We will present a further study of the effects of the metalHcity on 2MASS J05591914-1404488 
using a new set of synthetic models (del Burgo et al., in preparation). The new models will in- 
clude a completely new equation of state and (where possible and available) improved line data. 
Preliminary tests with the new equation of state show great improvements in the M dwarf regime. In 
addition, a physical model for the dust cloud formation coupled to the structure of the atmosphere 
is being developed; this is required at least for the transition T— >L— >M. 



5.2. Effective temperature 

For the two eailiest T dwarfs, 2MASS J05591914-1404488 and especially SDSSp J125453.90- 
012247.4, we found that the models used here cannot reproduce the observed spectra so well. The 



20Please give a shorter version with: \authorrunning and/or \titilerurming prior to \maketitle 



J 


1.0 
0.8 
0.6 








M=63 -j 


o 


0.* 










V 


0.2 
















1.205 


1.210 
wovelength [microns] 


1.215 


relotive flux 


1.0 
0.8 
0.6 
OA 
0.2 
















1.225 


1.230 
wovelength [microns] 


1.235 


relotive flux 


j:| 

o!9 

8-9 
















1.245 


1.250 
wavelength [microns] 


1.255 


relolive flux 


1.6 

1.0 
0.8 


u 






M = 60 "1 








1.265 


1.270 
wovelength [microns] 


1.275 


relotive flux 


1.6 
1.4 
1.2 
1.0 
0.8 

81 
















1.285 


1.290 

wovelength [microns] 


1.295 1.300 


relotive flux 


1.4 

i!o 

0.8 
0.6 
0.4 

6i 








M = 58 




1.J05 


1.310 


1.315 


1.320 



wovelength [microns] 



Fig. 11. 2MASS J15530228+1532369. Black and red lines correspond to the observed and mod- 
elled spectra, respectively. The model is for an effective temperature of 941 K and a surface gravity 
of lO^-* cm s-2. 

T4.5 dwarf 2MASS J05591914-1404488 has historically been considered enigmatic (see Vrba et 
al. 2004). Our best modelled spectra with < Teff > of 1002+278 K is, taking into account the 
errors, consistent with the estimate of 1231 K from the relation of Looper et al. and the value 
of 1200 K found by Gushing et al. (2008). Our best modelled spectrum of SDSSp J125453.90- 
012247.4 indicates an effective temperature above 2000 K, which is much higher than that found by 
Vrba et al. (2004) and Gushing et al. (2008), but consistent (i.e. within the errors) with a temperature 
around 1500 K. This value is similar to the estimate of 1370 K in the relation of Looper et al. 

We also find some differences for the latest T dwarfs. Our values of < Teff > for GL 570 D 
and 2MASS J0415 1954-0935066 are about 150 K higher than those from the literature (seefTTTi. 
Geballe et al. (2001) use R=400 spectra and accurate photometry of GL 570D to determine that 
Tetf of 784-824 K and log g=5.00-5.27 (cms^^), assuming an age of 2-5 Gyr. 

Fig. [16] shows Teff versus spectral type for the eight latest T dwarfs of our sample, and the 
Teff-spectral type relation found in Looper et al. (2008). Our values appear to be rather flat from 
T4.5 down to T8, which clearly contrasts with the trend delineated by Looper et al. (2008). We 
note, however, that such a trend is still compatible with our measurements if error bars are taken 
into account. The apparent constant temperature derived in our work may be due to a degeneration 
in the method (synthetic and observed dataset) that, despite the high spectral resolution of the 
observations, is not accurately sensitive to T^s and log g, partly due to the small wavelength range 
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Fig. 12. 2MASS J1217 1 1 10-03 11131. Black and red lines correspond to the observed and modelled 
spectra, respectively. The model is for an effective temperature of 922 K and a surface gravity of 
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coverage of the data. Gushing et al. (2008) found that the values of TeB obtained from J-band low- 
resolution spectra are consistent with those derived from fitting the full SED, using a model grid 
with temperature steps of 100 K. 

Fig. \T7\ leff) shows < T^ff > as a function of the colour J - [4.5] for the T dwarfs later than 
T4.5 (photometric data are compiled from Patten et al. 2006). These objects display a wide range 
of colours (~2 mag), however, our derived temperatures differ by less than 100 K. It seems difficult 
to reconcile such a large colour range with just one value of temperature, unless other atmospheric 
parameters (hke metallicity, cloud coverage, and others) are taken into account. 

5.3. Surface gravity 

Our values of log g are in good agreement with those of Knapp et al. (2004) (see Table fTTTl. whose 
logg values were derived from the comparison of observed H-K colors to model predicted H-K 
colors. The largest difference is for the T8 dwarf 2MASS J0415 1954-0935066, where the estimate 
of Knapp et al. (2004) is 0.7 dex higher than our value. Fig. [Vll right) shows log g versus infrared 
colour J - [4.5], where no obvious trend is apparent. 



22Please give a shorter version with: \authorrunning and/or \titilerurming prior to \maketitle 




1.2 JO 
wovelength [microns] 




1.250 
wavelength [microns] 




1.270 
wovelength [microns] 




1.J15 
wovelength [microns] 



Fig. 13. GL 570 D. Black line corresponds to the observed spectrum, and red and blue lines corre- 
spond to modelled spectra with the same ^=10^*^ cm s"^ and reif=948 K and 900 K. 



5.4. Mass and age 

The estimated masses in the T dwarfs of our sample are in good agreement with the model fit values 
of Burgasser et al. (2004) for three of the common objects (SDSS J162414.37+002915.6, 2MASS 
J15530228+1532369, 2MASS J0415 1954-0935066). The only exception is GL570 D. Burgasser 
et al. (2004) find a mass of <0.001 Mq, i.e., <1 Mj, and also provide a expected value of 30-50 
Mj. We find a value of 15 Mj, with an upper limit of 30 Mj that is the same as the lower expected 
value given by Burgasser et al. Our lower limit (6 My) is several times higher than the fitted value 
of Burgasser et al. (2004). 

Our result on the apparent young age of the field T dwarfs (see Table [TTI) is in agreement with 
recent kinematical studies based on proper motions and space velocities by Bannister & Jameson 
(2007) and Zapatero Osorio (2007). The latter authors found that about 40% of the L and T-type 
population of the solar neighborhood may have an age close to that of the Hyades cluster (around 
600 Myr), and that the brown dwarf population is kinematically younger than solar-type to early-M 
stars with likely ages in the interval 0.5-4 Gyr We also note that our upper limit of 2 Gyr for the 
widely studied object GL 570 D agrees with the lower limit assumed by Geballe et al. (2001). 
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Fig. 14. 2MASS J0415 1954-0935066. Black and red lines correspond to the observed and modelled 
spectra, respectively. The model is for an effective temperature of 947 K and a surface gravity of 
lO'*-^ cm s-^. 



6. Conclusions 

We conclude that the high resolution spectra corresponding to T dwarfs with spectral types 
later than T5 are well reproduced by the AMES-COND solar metallicity models provided by the 
PHOENIX code. The models reproduce in detail many faint absorption features in the high resolu- 
tion /-band spectra, which are mainly due to water vapor There are also strong K i lines, which turn 
out to be more difficult to model due to uncertain damping constants. The temperature and surface 
gravity determined from the comparison of the modelled and observed spectra are consistent with 
those found in the literature. We find a marked flat behaviour of T^ff with spectral type, although 
a possible gradient is compatible with the errors. High resolution spectroscopy {R ~ 20, 000) in 
the 7-band and AMES-COND models seems to be insufficient to show the existence of a possible 
gradient in effective temperature from early to late T dwarfs. The comparison between the spec- 
troscopically derived Te/f and logg of our targets and the evolutionary models by the Lyon and 
Arizona groups yields ages in the range 0.1-5 Gyr and masses between ~5 and 75 My for the target 
sample with spectral types >T5. For the earlier type dwarfs the spectral models do not provide 
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Fig. 15. Diagram T^g versus log g with the values we have obtained for our sample of T dwarfs 
(SDSSp J125453.90-012247.4 is not included) and the solar metalhcity models by Lyon (top) and 
Arizona {bottom). The models provide isochrones from 10 Myr to 10 Gyr (dashed lines) and curves 
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suitable fits, which is hkely due to the presence of condensate clouds that are not incorporated in 
the models. 
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